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ABSTRACT 

The importance of the radiative feedback from massive black holes at the 
centers of elliptical galaxies is not in doubt, given the well established relations 
among electromagnetic output, black hole mass and galaxy optical luminosity. 
In addition, feedback due to mechanical and thermal deposition of energy from 
jets and winds emitted by the accretion disk around the central black hole is also 
expected to occur and has been included in the work of several investigators. In 
this paper we improve and extend the accretion and feedback physics explored in 
our previous papers to include also a physically motivated model of mechanical 
feedback, in addition to radiative effects. In particular, we study the evolution of 
an isolated elliptical galaxy with the aid of a high-resolution 1-D hydrodynam- 
ical code, where the cooling and heating functions include photoionization and 
Compton effects, and restricting to models which include only radiative or only 
mechanical feedback (in the form of nuclear winds). We confirm that for Ed- 
dington ratios above 0.01 both the accretion and radiative output are forced by 
feedback effects to be in burst mode, so that strong intermittencies are expected 
at early times, while at low redshift the explored models are characterized by 
smooth, very sub-Eddington mass accretion rates punctuated by rare outbursts. 
However, the explored models always fail some observational tests. If we assume 
the high mechanical efficiency of 10~ 2 ' 3 adopted by some investigators, we find 
that most of the gas is ejected from the galaxy, the resulting X-ray luminosity is 
far less than is typically observed and little SMBH growth occurs. But models 
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with low enough mechanical efficiency to accomodate satisfactory SMBH growth 
tend to allow too strong cooling flows and leave galaxies at z = with E+A 
spectra more frequently than is observed. In a surprising conclusion we find that 
both types of feedback are required. Radiative heating over the inner few kpc is 
needed to prevent calamitous cooling flows, and mechanical feedback from AGN 
winds, which affects primarily the inner few hundred pc, is needed to moderate 
the luminosity and growth of the central SMBH. Models with combined feedback 
are explored in a forthcoming paper. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: 
active — galaxies: nuclei — quasars: general — galaxies: starburst 

1. Introduction 

All massive bulges and elliptical galaxies contain massive black holes at their center 
(hereafter SMBHs, e.g., see Kormendy & Richstone 1995, de Zeeuw 2001, Ferrarese & Ford 
2005), and when gas is added to the central regions for any reason the SMBH will accrete 
and emit energy. It is also clear that SMBHs have played an important role in the processes 
of galaxy formation and evolution (e.g., see Silk & Rees 1998; Fabian 1999; Burkert & 
Silk 2001; Cavaliere & Vittorini 2002; King 2003; Wyithe & Loeb 2003; Haiman, Ciotti & 
Ostriker 2004; Granato et al. 2004; Sazonov et al. 2005; Murray, Quataert & Thompson 
2005; Di Matteo, Springel & Hernquist 2005; Begelman & Nath 2005; Hopkins et al. 2006; 
Croton et al. 2006; Pipino, Silk & Matteucci 2008), as strongly supported by the remarkable 
correlations found between host galaxy properties and the masses of their SMBHs (e.g., see 
Magorrian et al. 1998, Ferrarese & Merritt 2000, Gebhardt et al. 2000, Yu & Tremaine 2002, 
McLure & Dunlop 2002, Graham et al. 2003, Marconi & Hunt 2003, see also Somerville 2008, 
Ciotti 2009). 

This basic fact leads immediately to a set of interrelated questions which must be 
addressed before we can understand either the resultant masses of the central SMBHs or 
the co-dependency of AGN feedback and galactic evolution, i.e. a) what physical processes 
control the SMBH accretion rate? For a given accretion rate, what form does the energy 
feedback take (photons, nuclear winds and jets)? b) How do these three forms of feedback, 
in turn, affect the accretion rate? 

In our prior papers (Ciotti & Ostriker 1997,2001,2007, hereafter C097, CO01 and CO07, 
and Ostriker & Ciotti 2005, hereafter OC05), we have focussed on the radiative, electromag- 
netic (EM) component, since that is most easily observed and our knowledge of it is most 
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certain (e.g., see Sazonov, Ostriker & Sunyaev 2004; Sazonov et al. 2007, 2008), and we al- 
ready provided indications that accretion feedback is the obvious solution to several related 
problems, such as the long lasting "cooling-flow" problem, and the fate of the large amounts 
of gas injected into the galaxy by the passively evolving stellar populations (e.g., Peterson 
& Fabian 2006, see also Binney 2001). This conclusion is also supported by other investi- 
gations implementing physically motivated feedback mechanisms (e.g., see Tabor & Binney 
1993, Binney & Tabor 1995, Omma et al. 2004, Churazov et al. 2005, Antonuccio Delogu 
& Silk 2008), and the computed solutions are characterized by relaxation oscillations (e.g. 
Cowie, Ostriker & Stark 1978; Milosavljevic et al. 2008). The general emerging picture is 
that energy output (radiative or mechanical) from the central SMBH pushes matter out, the 
accretion rate drops precipitously and the expanding matter drives shocks into the galactic 
gas. Then the resulting hot bubble ultimately cools radiatively (it is thermally unstable) and 
the consequent infall leads to renewed accretion; the cycle repeats, with the galaxy being 
seen alternately as an AGN/starburst for a small fraction of the time and as a "normal" 
elliptical hosting an incipient cooling catastrophe for much longer intervals. Nowadays, sev- 
eral observations support the finding that accretion on central SMBHs is in fact a highly 
unsteady phenomenon (e.g., see Martini 2004, Goncalves, Steidel & Pettini 2007; Prochaska 
& Hennawi 2008, Hopkins & Hernquist 2008). 

However, although our previous work has stressed the significance of radiative heating 
near to the SMBH (in addition to the well known Eddington momentum input, e.g. see 
Dorodnitsyn, Kallman, & Proga 2008; Shi & Krolik 2008), other investigators (e.g., Binney 
& Tabor 1995, Tabor & Binney 1993, Begelman & Nath 2005, Begelman & Ruszkowski 
2005, Di Matteo, Springel & Hernquist 2005, Springel, Di Matteo & Hernquist 2007) have 
focussed on the also important but highly uncertain mechanical feedbaclo An obvious 
source of mechanical feedback is certainly represented by the radiatively driven winds from 
the broadline regions (BLRs), whose parameters are well observed (e.g., Chartas et al. 2003, 
2007; Crenshaw, Kraemer, & George 2003; Blustin et al. 2007; Hamann et al. 2008), so 
that their energy, momentum, and mass input can be added to the codes in a fairly direct 
way. Although the overall efficiency of such inputs appears to be modest, this component of 
feedback couples with great effectiveness to the ambient gas. In addition, highly collimated 
jets, especially at low Eddington ratios are observed to put out energy in amounts comparable 
to the EM output, but it is not clear how efficiently such narrow jets can couple to the ambient 
fluid (see Vernaleo & Reynolds 2006). 



1 Note that also a purely radiative feedback produces ultimately some form of mechanical feedback, in the 
form of shock waves, e.g. CO01 and CO07, or in the form or radiation driven winds as found in 2D and 3D 
simulations (Proga 2007, Proga et al. 2008b, Kurosawa & Proga 2008ab). 
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The purpose of the current paper is to introduce a physically based modelling of mechan- 
ical energy input to our code, and to supplement the detailed treatment of radiative effects, 
in a fashion that is consistent with current theory and guided by observations. At the same 
time, we further improve the galaxy models used for the simulations, both in the stellar and 
dark matter distributions, following the results of the latest observational works. In par- 
ticular, in this paper we restrict our exploration to purely radiative and purely mechanical 
feedback, in order to better understand the specific properties of the two mechanisms when 
considered separately. Other aspects of purely mechanical models are investigated in detail 
in Shin, Ostriker & Ciotti (in preparation, hereafter Paper II), while we reserve to a third pa- 
per (Ciotti & Ostriker, in preparation, hereafter Paper III) the discussion of models in which 
both radiative and mechanical feedback effects are included in the hydrodynamical code. 
Finally, the X-ray observational properties of the most successful models will be presented 
in Pellegrini, Ciotti & Ostriker (in preparation, hereafter Paper IV. See also Pellegrini, Ciotti 
& Ostriker 2009), to be compared with the results of observational works (e.g., O'Sullivan, 
Ponman, & Collins, 2003; Diehl & Statler 2008). 

We continue to find that the situation is inherently complex. Mechanical energy output 
of the AGN is lower than the EM output but couples more effectively to ambient gas and 
thus is optimal in the innermost few hundred parsecs in shielding the SMBH from excessive 
accretion. But the radiation heating input to the galaxy is carried primarily by long mean 
free path (X-ray) photons and is most effective in heating the inner several kpc and thus 
moderating or forestalling cooling catastrophes in the ambient ISM of the galaxy. We tenta- 
tively conclude the BLR winds do have an essential effect in regulating SMBH accretion and, 
with the radiation heating feedback, cooperate in removing significant amounts of gas from 
giant ellipticals (see Schawinski et al. 2008) but that the jets largely escape from isolated 
galaxies. However, for Brightest Cluster Galaxies, the energetic jets can provide significant 
input to the cluster gas - thereby helping to prevent or moderate cluster cooling flow catas- 
trophes (e.g., Voit & Donahue 2005; Peterson & Fabian 2006; Rafferty, McNamara & Nulsen 
2008; McCarthy et al. 2008). 

The paper is organized as follows. In Section 2 we describe how the new galaxy models 
adopted in the simulations are built, the details of the input physics (focusing in particu- 
lar on the physics of accretion disk and on the mechanical feedback), and their numerical 
implementation. In Section 3 we present two different classes of models, both of them pre- 
senting some interesting aspects but overall failing to reproduce the situation observed in 
real galaxies. We thus have a first class of models in which only radiative heating is taken 
into account, and a second class in which only mechanical feedback is considered. Finally, 
in Section 4 we discuss the main results obtained. 
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2. The models 

The galaxy models and the input physics adopted for the simulations have been improved 
with respect to our previous explorations. In the following we describe the new ingredients 
of the current models, while for the unchanged parts we refer to our previous papers: a 
comparative summary of the present and past treatments is given in Table 1. 



2.1. Structure and internal dynamics 

In C097 and CO01 the galaxy models utilized a King (1972) stellar distribution plus a 
quasi-isothermal dark matter halo, in line with the models then used for cooling-flow studies. 
However, the existence of large cores of constant surface brightness has been clearly ruled 
out, as HST observations have shown how the central surface brightness profile is described 
by a power-law as far in as it can be observed, i.e. to ~ 10 pc from the center for Virgo 
ellipticals (e.g., see Jaffe et al. 1994, Faber et al. 1997, Lauer et al. 2005). In Pellegrini & 
Ciotti (1998) and in CO07 a stellar density distribution described by the more appropriate 
Hernquist (1990) model has been adopted, and also the dark matter halo was described by 
an Hernquist profile, which is quite similar in its central regions to the dark matter halos 
obtained from cosmological simulations (e.g., Dubinski & Carlberg 1991; Navarro, Frenk & 
White 1997; Fukushige & Makino 1997). 

However, the Hernquist profile is characterized (as the de Vaucouleurs i? 1 / 4 density 
profile) by a sizable central depression in the isotropic velocity dispersion, which is not 
observed. This problem is now fixed. In fact, the velocity dispersion of the Jaffe (1983) 
profile is monotonically decreasing. In addition, observations support the idea that ellipticals 
are characterized by a total density distribution well described by a r~ 2 profile (e.g., see Treu 
& Koopmans 2002,2004; Rusin et al. 2003, Rusin & Kochanek 2005, Koopmans et al. 2006, 
Gavazzi et al. 2007, Czoske et al. 2008; Dye et al. 2008). For these reasons, here we adopt 
Jaffe stellar models embedded in a dark halo so that the total mass profile decreases as r~ 2 . 
Therefore, the stellar density profile of the galaxy models is 

M*r* m 
P *~47rr 2 (r, + r) 2 ' 1 j 

where M* and r* are the total stellar mass and the scale-length of the galaxy, respectively; 
we recall that for the Jaffe model the effective radius is R c = r*/0.7447. The total density 
profile is then given by 



Table 1. Synoptic table of feedback model development 



Input Physics 


C097 


CO01 


OC05 


CO07 


Papcrl [Radiative] 


Paperl [McchamcalJ 


Paperlll 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Galaxy model 


Kimr+Q T H 




Kine+O T H 


TT p r ti n 1 1 i ^ t -\- 1~T p v ti n 1 1 i ^ t 


laffe in S T S 


Taffe in S T S 

O Cill^ 111 k-J ♦ _L . Lj ♦ 


laffe in S T S 


Star formation 


No 


No 


No 


Yes 


Yes 


Yes 


Yes 


Circumnuclear disk 


No 


No 


No 


Yes 


Yes 


Yes 


Yes 


ADAF 


No 


a Yes 


No 


No 


Yes 


Yes 


Yes 


Heating and Cooling 
















Compton 


b Yes 


b Yes 


Yes 


Yes 


Yes 


c Only cooling 


Yes 


Photoionization 


No 


No 


Yes 


Yes 


Yes 


c Only cooling 


Yes 


Radiation Pressure 
















e~ scattering 


Yes 


Yes 


Yes 


Yes 


Yes 


d No 


Yes 


Photoionization 


No 


No 


No 


Yes 


Yes 


d No 


Yes 


Dust 


No 


No 


No 


Yes 


Yes 


d No 


Yes 


Mechanical Feedback 
















AGN wind 


No 


No 


No 


d No 


d No 


Yes 


Yes 


AGN jet 


No 


No 


No 


No 


d No 


d No 


d No 



Note. — Different names indicate the papers as in the text. Paper I indicates this paper; Paper III; Q.I.H.= Quasi-Isothermal Halo; S.I.S.= Singular Isothermal 
Sphere: in the present models the stellar density is immersed in a dark matter halo so that the total density profile is proportional to r~ 2 . 

a Only a few test models were computed. 

b In C097 a very high Compton temperature Tx for the emitted accretion luminosity spectral distribution was used, while in CO01 a large range of values for 
Tx was explored. Starting from OC05 the value of Tx ~ 2 10 7 K was fixed according to the observational estimates of Sazonov, Ostrikcr & Sunyaev (2004). 

c Heating computed but not added to the hydrodynamics. 

d Computed but not added to the hydrodynamics. 
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where the dimensionless factor 1Z > 1 controls the amount of dark matter contained within 
the half-mass radius of the stellar component. With this choice, the circular velocity of the 
model is given by v c = a/ G1ZM*/r*. The presence of significant amounts of dark matter 
in the central regions of normal galaxies remains controversial (e.g., see Binney & Evans 
2001), but in general observational studies in the optical (e.g., Saglia et al. 1993; Bertin 
et al. 1994; Cappellari et al. 2006, Douglas et al. 2007) and in X-rays (e.g. Fabian et al. 
1986, Humphrey et al. 2006) seems to indicate that dark matter begins to be dynamically 
important at 2 — 3R e . For these reasons we fix 1Z to the minimum admissible value 1Z = 1, 
corresponding to the same amount of dark and visible matter within the spatial half-mass 
radius of the stellar distribution. 

All the dynamical and phase-space properties of the resulting two-component galaxy 
models are given in Ciotti, Morganti & de Zeeuw (2008), and here we report only the 
formulae for the quantities of interest: in particular, the model central projected velocity 
dispersion a Q (obtained by solving and projecting the Jeans equations under the assumption 
of orbital isotropy) is given by 



0"o = 



y/2' 



(3) 



Note that we do not consider the effect of M BH on a , in accordance with estimated values 
for the radius of the SMBH sphere of influence (e.g., Riciputi et al. 2005). The parameters 
describing the galaxy model are determined following CO07, i.e., first we assign a value for 
the central velocity dispersion a , and we determine the galaxy present day blue luminosity 
Lb and effective radius R c from the Faber-Jackson and the Fundamental Plane relations 
(eqs. [3]-[4] in CO07). We then fix the relative amount of dark matter to stars within R e by 
assigning 71, thus determining the total stellar mass of the galaxy and finally the stellar 
mass-to-light ratio Y* = M*/Lb- 

An important ingredient in the energetics of the gas flows, namely the thermalization 
of the stellar mass losses due to the stellar velocity dispersion (e.g., see Parriott & Bregman 
2008), depends on the radial trend of this latter quantity which, for the isotropic models 
here considered is given by 



p*cr* = p*a t0 + 



GM.Mbh 



1 -2s + 6s 



12s A 1 

— 4 In 1 + - 

3s 3 1 + s) \ s 



s = 



(4) 



where a* Q is the isotropic 1-dimensional stellar velocity dispersion without the contribution 
of the SMBH (Ciotti et al. 2008). Note that, at variance with the estimate of a , in the 
thermalization of the stellar mass losses we also consider the (time increasing) contribution 
of the gravitational field of the central SMBH. 
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2.2. The circumnuclear disk and the SMBH accretion luminosity 

At the onset of the cooling catastrophe a large amount of gas suddenly flows onto the 
central regions of the galaxy, and this induces star formation and accretion on the central 
SMBH, producing a burst of energy from the galaxy center. However observations of our 
own galactic center and high resolution studies of other nearby systems indicate that, in 
addition to the central starburst with radius ~ 100 — 300 pc, accretion onto the central 
SMBH is mediated by a small central gaseous disk within which additional significant star 
formation occurs, and the remaining fraction of gas either is blown out in a BLR wind, or it 
is accreted onto the central SMBH. In our treatment the disk is not simulated with hydrody- 
namical equations, but its description is needed to obtain important quantities required by 
the code. In an improvement over CO07, we now also model and add to the hydrodynamics 
the mechanical feedback produced by the disk wind (see Table 1), but for future reference, 
in the following we also describe the modeling of a nuclear jet. However, while the jet con- 
tribution to the circumnuclear disk mass balance is considered in the numerical integration, 
the associated mechanical feedback is not added to the current hydrodynamical simulations. 

The circumnuclear disk, which is the repository of the gas inflowing at a rate Mf s 
from the first active mesh point R\ of the hydrodynamical grid, and which feeds the central 
SMBH at a rate Mbh, contains at any time the mass gas M dg and a total stellar mass 
M d * = M d i* + Mdh*, which is divided among low and high mass stars (with the division mass 
at 8M Q ). The disk also contains a mass M rem of remnants from the earlier generations of 
evolved stars. 

In the adopted scheme the accretion rate on the central SMBH is given by 

M BH = — , (5) 



where 



M dg _ Mm • Z/Edd 



r d 2M Edd e oc 2 

are the fiducial depletion rate of gas from the circumnuclear disk, its normalized value, 
and the Eddington mass accretion rate, respectively. The reference radiative efficiency €q is 
defined in eq. (IE])- Equations ([5])- ([6]) are designed to guarantee that when rjd <C 1 the gas is 
accreted onto the central SMBH at the rate M fid , while Mbh = 2M Edd for r/ d ^> 1 (i.e., we 
allow for possible moderate super-Eddington accretion; note however that outside the first 
grid point R± the flow accretion rate is limited in a self-consistent way by feedback effects). 
From eq. ([5]) we calculate the instantaneous bolometric accretion luminosity as 



e EM M BH c 2 , (7) 
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and, at variance with CO07 (see also Table 1), here we adopt an "ADAF-like" radiative 
efficiency 

Am . M B h , s 

where A is a free parameter so that e^M ~ eoArh for m <C A" 1 . In our simulations we fix 
A = 100 (see, e.g. Narayan & Yi 1994, and CO01, where a very preliminary investigation 
of ADAF effects on radiative feedback was carried out), and we adopt eo = 0.1 or 0.2 (e.g., 
see Noble, Krolik & Hawley 2008). As usual in accretion theory, we finally introduce the 
normalized accretion luminosity 

_ Lbh = Am 2 
L Edd 1 + Am 

where the last expression derives from the ADAF phenomenological description. 

There are a few lag times in our problem which are expressed as follows. The first is the 
instantaneous disk lag time, appearing in eq. ([6]), and already considered in CO07, defined 

as 

where a ~ 10~ 2 — 10 _1 is the disk viscosity coefficient, and and Mbh are the instantaneous 
values of the fiducial radius of the circumnuclear disk and the mass of the central SMBH. 
In CO07, the disk radius was maintained fixed to Ri, instead now we use the scaling 
predicted by thin-disk theory 

R d (t) = f d R 1 x (_*L , (11) 

(e.g., see Morgan et al. 2007) where Mbho is the central SMBH mass at the beginning of 
the simulation. We assume fd = 0.4, so that i?d(0) — 2 pc for an initial SMBH mass of 
~ 1O 8 M . 

The second characteristic time (that was not considered in CO07) is the instantaneous 
infall lag time from Ri to the disk, estimated as 



Rx 2GM BH 

^ = — , v s = W — , 12 

v s V Ri 

so that the effective rate at which gas accretes on the disk is obtained by solving the differ- 
ential equation 

dMf M x - Mf , . 
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where M\ is the instantaneous rate at which gas flows through the first active grid pointEl 
It follows that when M\ provided by hydrodynamics drops to zero the circumnuclear disk 
experiences a fueling declining exponentially with time. 

The disk total gas mass M dg is not only the source of SMBH accretion, but also of star 
formation in the disk: we assume that a fraction of M dg is converted into stars at a rate 
77*Mfid (where r]* ~ lOM^/M^n), and that another fraction of M dg is lost as a disk wind and 
as a jet at instantaneous rates given by ^ w M B h and rjjM B n, so that the equation for the gas 
mass in the disk is 

^5 = Mf - (1 + Vw + t^Mbh - V*M m . (14) 

Note that the presence of a jet was not considered in CO07, so that r/j = therein. The stars 
formed in the disk are described separately as a function of their mass, i.e., high- mass stars 
(M > Mn = 8M ) produce a total disk mass Mdh*, and low-mass stars (M in f < M < Mn) 
contribute to a disk mass Mai* according to the equations 

dM dU . M dU dM dh * ■ M dh * 

— 37— = (i - fh)v* M &d ; —77— = Uv* M &d • (15) 

at r*i at r* h 

For the characteristic evolutionary times we adopt r*i = r opt and r*h = rji given in CO07, 
while we assume /h = 0.5, corresponding to a top-heavy Salpeter-like initial mass function 
(e.g., see Nayakshin & Sunyaev 2005, Nayakshin et al. 2006) of slope x ~ 1.16 and minimum 
mass Mi n f = 0.1M o . The associated optical (£d,o P t) and UV (£d,uv) luminosities of the stellar 
disk are calculated following the scheme described in CO07. Finally stellar remnants mass 
in the disk evolves as 

dM rem _ M dU Mdh, rifi , 

j, Jrem,l i Jrcm,h ) \ 

etc r*i r^h 

where / rem ,i = 0.2, / rem) h = 0.09. The equation for the mass loss associated with the disk 
wind is 

dM dw ■ M d i* Md^ 

— — = T] w M B n + (1 - / rem ,l) h (1 ~ /rem.h) : (17) 

at r*i r*h 

the first term is a mass loss driven as a wind by the central SMBH, and the second and third 
are from high mass and low mass stars in the central disk. 

At variance with CO07, where the effects of mechanical feedback were not considered 
and the wind treatment was extremely simplified, here we explore two different classes of 
models, that we call Type A and Type B. In particular, we use for the ratio of the wind 



2 Mi is taken positive in case of accretion and zero in case of outflow at Ri. 
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Fig. 1. — Mechanical feedback properties as a function of Eddington ratio / with ADAF 
coefficient A = 100 (eq. [9]). Red solid lines refer to wind properties of Type B models with 
maximum wind efficiency e^f = 5 10~ 3 and 77^ = 1800e^ (eq. [23]). From bottom to top pan- 
els: ejected mass fraction (eqs. [IB]- [23]); wind mechanical efficiency (eq. [2"T]); wind velocity 
(eq. [22]); normalized wind momentum per unit solid angle P/Pedd = \/2e w ?7 w m/ (2Af2 w ) 
(where PEdd = M^ddC, see eq. [26]). The wind opening angle is given in eq. (I3TT) . Black solid 
lines refer to wind properties of Type A models, where e w = 1.5 10~ 4 , 7/ w = 2, AQ W = ir. 
Dotted lines represent the corresponding jet properties, parametrized as in eqs. (fl9l) and (1241) - 
(I26p . with AQ W = 2.5 10 -2 . Finally, the green line is the ADAF modulated EM efficiency as 

civpn hv pn flSll with f n = fl 1 anH A = 1DD 
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outflow rate to the SMBH accretion rate 




4 1 + 0.25/' 



[A] 
[B]. 



(18) 



Therefore, Type A models correspond to the models in CO07. In case of Type B models, 
whose parameterization is introduced to mimic the results of 2D simulations, the efficiency 
of ejecting a wind increases with increasing Eddington ratio, rj^ fixes the maximum value 
of rj w , and the factor 3/4 takes into account the maximum possible value for the scaled 
accretion luminosity /. The situation is illustrated in the bottom panel of Fig. 1, where the 
solid lines represent the coefficient i] w for Type A models (black) and Type B models (red, 
with 77^ = 1800e^). As a consequence of this choice, winds in Type A models eject more 
mass at low accretion luminosities, while the opposite happens in Type B models for />l/3. 

In particular, our scaling of the radiation driven wind mass loss rate and mechanical 
power are motivated both by current simulations and by observations of massive outflows 
from quasars (e.g., Ganguly & Brotherton 2007; Holt, Tadhunter & Morganti 2008). Some 
quasars show broad absorption lines (BALs) which are the most dramatic evidence for winds 
in AGN. BALs are almost always blueshifted relative to the emission-line rest frame, indi- 
cating the presence of outflows from the active nucleus, with velocities as large as 0.2 c (e.g., 
Turnshek 1998). BALs are observed not only in the UV but also in the X-rays: for example, 
Chartas, Brandt & Gallagher (2003) discovered a very broad absorption line in the X-ray 
spectrum of PG 1115+80. It is commonly accepted that these mass outflows in quasars are 
very likely disk winds driven by radiation from accretion disks (e.g., see Konigl 2006, Proga 



Of course, determining the mass loss rate and mechanical power based on observations 
requires modelling, because observations do not directly provide the information about the 
wind column density and photoionization structure (e.g., Arav et al. 2007, and references 
therein). In addition, one must assume something about the wind covering factor. Therefore, 
a physical model of disk winds is needed to complement data analysis and to estimate the 
key wind properties. In particular, hydrodynamical simulation of radiation driven disk winds 
allows us to explore the impact upon the mass-loss rate, Md w and outflow geometry caused 
by varying the system luminosity and the radiation field geometry. For example, Proga at 
al. (1998) showed that winds driven from, and illuminated solely by, an accretion disk yield 
complex, unsteady outflow. In this case, time-independent quantities can be determined only 
after averaging over several flow timescales. On the other hand, if winds are illuminated by 
radiation mainly from the central engine, then the disk yields steady outflow. Proga et al. 
(1998, see also Proga 1999) found that Mdw is a strong function of the total luminosity, while 
the outflow geometry is determined by the geometry of the radiation field. In particular, 
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for a relatively small I, there is no wind but only a puffed up disk. Then as I increases a 
strong equatorial wind develops. For very high I, the wind becomes bipolar. For high system 
luminosities Mdw of the disk scales with the luminosity in a way similar to stellar mass loss 
rate, i.e., Mdw oc / L7 (see Castor, Abbott & Klein 1975). As the system luminosity drops 
below a critical value [~ 2LEdd/(l + M max )], M& w decreases quickly to zero (M max is the total 
line opacity for an optically thin case). 

The wind power and geometry depend on the luminosity and also on the mass of SMBH. 
Proga & Kallman (2004, see also Proga et al. 2000) showed that for M BH = 1O 8 M and 
/ ~ 0.5 a strong wind develops whereas for I ~ 0.1 there is no disk wind (eq. [18] for Type 
B models reflects this trend). The primary reason for this luminosity sensitivity is the fact 
that the mass flux density of the wind decreases strongly with decreasing disk luminosity 
and the wind is more subject to overionization. Proga & Kallman (2004) also found that for 
a fixed Eddington fraction (e.g., 50%) it is easier to produce a wind for Mbh > 1O 7 M than 
for Mbh < 1O 7 M . This result is a consequence of the decrease of the UV contribution to 
the disk total radiation with decreasing mass of the BH in the Shakura & Sunyaev (1973) 
disk model for a fixed Eddington fraction. 

In addition to the treatment of CO07, we now consider another mass component ejected 
by disk, i.e. a nuclear jet with instantaneous mass flow 

dM \ »v 0.2 

— ^ Mbh ' ^ = (ITTbW ; (19) 

the trend of rjj as a function of / is shown in the bottom panel of Fig. 1 with the dotted line: 
it is apparent how the mass ejected by the jet is always negligible with respect to the wind 
mass loss in Type A models, while it is slightly dominant over the wind in Type B models 
at low luminosity ratios. 

We stress that, even though we consider this additional mass component in the circum- 
nuclear wind mass balance (eq. [H]), the associated feedback effects on the galaxy ISM are 
not taken into account in the simulations (see Table 1), and for the moment we are assuming 
that the mass, momentum, and energy fluxes in the jet escape the galaxy. In the code, all 
the equations presented in this Section are integrated numerically with a first order finite 
difference scheme. 



2.3. The mechanical feedback treatment 



We now discuss how the kinetic energy, momentum and mass of the BLR wind are 
transferred to the ISM. For sake of completeness, we also describe the jet treatment. 
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As in CO07, the fiducial instantaneous mechanical luminosity of the disk wind is given 

by 

L dw = e w M BH C 2 + e n C 2 (l - / renll h)— (20) 

where e w is the mechanical efficiency of the wind, and the second term describes the energetic 
associated with the SNII explosions of the high-mass stars in the circumnuclear disk. In 
analogy with eq. (TTBl . we assume 



~W 5 

M 



e- [A] 



= < 3g j [B1 (21) 

4 1 + 0.25/' [ J ' 

In Type A models we explore the range 3 10~ 5 < e^ 1 < 5 10~ 3 (see Table 2), i.e also values 
significantly lower than adopted by Hernquist and collaborators (5 10~ 3 , see Di Matteo et al. 
2005). For reference, in CO07 (where e w was also kept constant), we adopted e^f = 5 10~ 4 , 
but the resulting mass, energy and momentum sources were not added to the code. In 
case of Type B models, where the wind efficiency is a function of the normalized accretion 
luminosity, is the maximum possible value (reached for 1 = 2), and the different values 
adopted in the simulations are listed in Table 2. In both cases the instantaneous disk wind 
velocity is given by 

~ f^ c (22) 
Mdw V ' 

where the last expression neglects the mass return contribution of massive stars in the cir- 
cumnuclear disk (see eq. [20]). Note that in Type A models v w is in the range 2 10 3 — 2 10 4 
km s -1 (as a function of the specific assumed value for e w ), in agreement with observations 
of BLRs (e.g., Crenshaw et al. 2003). For the same reasons, in Type B models we require 
v w = 10 4 km s _1 , so that rj^ and e^f are linked by the relation 

rg = 1800e^. (23) 

In analogy with the wind component, the instantaneous jet mechanical luminosity is 
written as Q125 

L^qAW, q=( ITW , (24) 



and the jet velocity is given by 




(25) 



which, for our chosen parameterization gives high but subrelativistic jet velocity of v-Jc ~ 
^q-i.65 f or [>q i ( see pjg_ x). Note that from eqs. (21) and (24) jets are far more efficient 



- 15 - 



than nuclear winds at low Eddington ratios (Ghisellini 2008, Jolley & Kuncic 2008). The 
dependence of wind and jet velocities on the normalized accretion luminosity is shown Fig. 1 
(second panel from the top), where solid lines refer to the wind and the dotted line to the 
jet. Finally, the wind and jet momentum are defined as 

m-j = MjVy, m w = M dw v w . (26) 

We now illustrate how we distribute the mechanical feedback over the galaxy ISM. In the 
following we describe the procedure both for the wind and for the jet, even though in the 
hydrodynamical equations we take into account only the wind component. First we introduce 
the instantaneous wind and jet lag times 

r wj ee ^ (27) 

from the center to the first active grid point R\ (where the subscript indicates the specific 
component - disk wind or nuclear jet - considered), and at each time step we compute the 
time-lagged values for mass, momentum, and kinetic energy at R± by solving the differential 
equation 

dX l X-X l 



dt 



(2? 



where X[ is the generic lagged variable associated with the instantaneous unlagged value X. 
Outside Ri we then distribute mass, momentun and kinetic energy over the hydrodynamical 
grid (outside Rx), by integrating numerically the phenomenological differential equation 

d\nY wj _ Pi SM (r) r d\nY wi 
d\nr Pwj(r) v wj dt 

where Y w j is the mass, momentum and energy of the disk wind/jet component at distance r 
from the center, -P w j( r ) is the local wind/jet pressure, and for each quantity Y(Ri) = X t . In 
this paper we restrict to simulations where the time derivative is neglected, while in Paper 
III also this latter term is considered. In practice, we first integrate eq. (j2"9~|) for the wind/jet 
pressure, i.e., 



where Y w j is the effective wind/jet momentum crossing the shell of radius r, so that eq. (1291) is 
a non-linear differential equation for Y W y Once the equation is integrated, the radial behavior 
of P W j and the r.h.s. of eq. ( 1291) are known over the whole grid, and the equation can be 
integrated for mass and energy. We thus determine how much of the mechanical energy, 
momentum and returned mass are deposited in the ISM at each radius. 
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A physical justification of eq. fl29l) is given in Appendix A, but it can be thought of 
phenomenologically in the following way. If the pressure corresponding to the momentum 
flow within the jet or wind is much greater than the pressure in the ambient gas, very little 
mass, momentum and kinetic energy is taken from it and deposited in that ambient gas. 
But when the r.h.s. of eq. (1291) approaches unity, the "working surface" has been reached 
and the jet or wind discharges its content. 

The solid angle in the denominator of eq. f)30p is the opening angle of the wind and of 
the jet, and the factor of 2 accounts for the biconical nature of the flow. While for the jet 
we assumed in all the simulations the fiducial value Af2j = 2.5 10~ 2 , for the wind case we 
adopt 

AQw= {lmin(VW+^,l), [B], (31) 

where case B is designed to mimic the behavior found in radiation driven winds: higher 
luminosity corresponds to a larger opening angles. The constant inside the square root is 
fixed to a = AQ^/tt, so that for small values of accretion luminosity the wind opening angle 
coincides with the jet opening angle. Finally, note that the almost linear dependence of AQ W 
on I for I > a assumes that the linear opening angle depends on \fl for this regime. The 
wind/jet momentum (with the opening solid angle factor) as a function of I is shown in Fig. 1 
(top panel), and it is apparent that this quantity is dominated by the wind contribution (solid 
lines) for Z>1(T 2 . 

To implement numerically the mechanical feedback terms, we finally compute the nu- 
clear wind mass, momentum and kinetic energy per unit volume deposited in each shell 
as 



Source ■ - 3 Y wj(Ri) Y vi (Rj +1 ) 
bource^-^ R s +i _ R 3 



and we add them (only for the wind component) to the r.h.s. of eqs. (56)- (58) in CO07. In 
practice, we assume that the quantities of mass, momentum and kinetic energy transferred 
to the ISM within the opening angle are finally distributed over the whole solid angle, a 
forced choice for our ID code. 

Summarizing, our simple formulae for Type B models capture some of the key wind 
properties of quasar winds: the wind opening angle, mass loss rate, and power increase with 
the accretion luminosity, and the fact that there are negligible radiation pressure driven 
winds for /<0.01. Figure 1 shows (going downwards from the top) our modeling of the 
inputs of momentum, velocity, energy and mass, respectively. The simplest case of constant 
wind mechanical efficiency (Type A) is shown as a solid black line, whereas the one in better 
accord with radiatively driven winds (Type B) is shown in red. A hypothetical relativistic jet 
model ( "radio mode" ) is given by the dotted lines and the green line indicates the ADAF-like 
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treatment of the radiative efficiency. For Eddington ratios above 0.01 the solution represents 
the typically observed (and computed) "AGN mode" of high radiative efficiency and low 
mechanical energy output efficiency, whereas at low Eddington ratios (Z^10 _L5 ) the output 
shifts to radio mode of low EM output and high mechanical energy efficiency in a relativistic 
jet. As we will see most of the SMBH growth occurs in the AGN mode and so the net (mass 
weighted) efficiencies are in agreement with well known estimates (e.g., Soltan 1982, Yu & 
Tremaine 2002, Merloni & Heinz 2008, Martinez-Sansigre & Taylor 2008), but much of the 
energy output occurs in the radio mode with important consequences for the energy balance 
in clusters of galaxies (e.g., see Mc Namara & Nulsen 2007, and references therein). 



2.4. The unchanged input physics 

The stellar mass loss rate and the SNIa rate associated with the initial stellar distribution 
are the main ingredients driving evolution of the models. In the code the stellar mass losses 
- the source of fuel for the activity of the SMBH - follow the detailed prescriptions of the 
stellar evolution theory, and we use exactly the same prescriptions as in CO07 (see Sects. 2.2 
and 2.3 therein). 

The radiative heating and cooling produced by the accretion luminosity are numeri- 
cally computed as in CO07 by using the Sazonov et al. (2005) formulae, which describe 
the net heating/cooling rate per unit volume E of a cosmic plasma in photoionization equi- 
librium with a radiation field characterized by the average quasar Spectral Energy Distri- 
bution derived by Sazonov et al. (2004, see also Sazonov et al. 2007, 2008), whose as- 
sociated spectral temperature is Tx — 2 keV. In particular, Compton heating and cooling, 
bremsstrahlung losses, line and continuum heating and cooling, are taken into account. Also 
the star formation over the galaxy body, the radiation pressure due to electron scattering, 
to phoionization, and finally to UV, optical and infrared photons on dust, are treated as in 
CO07, where the derivation and the numerical integration scheme of the radiative transport 
equations is described in detail. Finally, all the relevant information about the numerical 
code and the hydrodynamical equations can be found in CO01 and CO07. 



3. Model evolution 

We now show the main properties of representative models characterized by a stellar 
mass M* = 2.9 x 1O 11 M , a blue optical luminosity L B = 5 10 10 L B © (corresponding to a 
stellar mass-to-light ratio in the blue band of ~ 5.8), a Fundamental Plane effective radius 
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R e = 6.9 kpc, and a central aperture velocity dispersion cr = 260 km s" 1 . The stellar 
distribution is immersed in a dark matter halo so that the total mass density distribution is 
proportional to r™ 2 , and an identical amount of stellar and dark matter is contained within 
the stellar half-mass radius. The initial SMBH mass follows the present day Magorrian 
relation, with Mbh — 10 _3 M*, as it is believed that the bulk of the SMBH mass is assembled 
during the process of galaxy formation (e.g., Haiman, Ciotti & Ostriker 2004; Sazonov et al. 
2005), but this process is not modeled in the present simulations. Therefore, these models are 
not appropriate as initial conditions for cosmological simulations, because their parameters 
are fixed to reproduce early-type galaxies similar to those observed in the local universe 
(at z = 0), and also because we set outflow boundary conditions at the galaxy outskirts 
(~ 250 kpc): from this point of view, the simulations represent an isolated elliptical galaxy 
(consistently we are not considering the effects of possible merging on the galaxy evolution). 
A central cluster galaxy would have more difficulty generating winds and would suffer from 
bursts of cluster gas inflow. We adopted this framework to adhere to the standard approach 
followed in "cooling-flow" simulations. In future explorations we will address in a more 
consistent way the problem of the galaxy structural and dynamical modifications due to 
star formation and mass redistribution over an Hubble time, and the compatibility of the 
obtained galaxies with the present-day scaling laws of elliptical galaxies. We also stress that 
the models presented are just a representative sample out of several tens of runs that have 
been made, characterized by different choices of the input parameters (often outside the 
currently accepted ranges). 

The initial conditions for the ISM are represented by a very low density gas at the local 
thermalization temperature. The establishment of such high-temperature gas phase at early 
cosmological times is believed to be due to a "phase-transition" when, as a consequence 
of star formation, the gas-to-stars mass ratio was of the order of 10% and the combined 
effect of SNIa explosions and AGN feedback became effective in heating the gas and driving 
galactic winds (e.g., Sazonov et al. 2005). Several theoretical arguments and much empirical 
evidence, such as galaxy evolutionary models and the metal content of the Intra Cluster 
Medium support this scenario (e.g., Renzini et al. 1993; OC05; Di Matteo et al. 2005). For 
the reasons above, in the simulations here presented (as well as in all others simulations not 
shown), we assume that the galaxy stellar component at the beginning of the simulation 
is 2 Gyr old, and the simulations span 12 Gyr, so that the cosmic time at the end of the 
simulations is 14 Gyr. 

Important quantities associated with the model evolution are the mass (luminosity) 
accretion weighted EM and mechanical efficiencies, defined in natural way as 



< ^EM > = 



/ 6emM B h^ 



/ e w M B udt 
AM BH 



(33) 



AM BH 
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where AM B h is the SMBH accreted mass over the time interval considered. In case of Type 
B models we also compute the luminosity weighted average fraction of the sphere covered 
by the wind, defined as 

< AO, > s S^f-f , (34) 
4tt / L BK dt v 1 

At variance with eq. (33) we restrict the computation of the integrals to high luminosity 
phases only, that we fiducially assume defined by Lbh > O.lLEdd- Following CO07, we 
finally compute for the various luminosities the duty cycle over a period of time At 

_ UU, Ld(f 

- At S :_ Al L^ (35) 

with this definition a square wave with fraction of time / in the high state would have 

/duty = /• 

We summarize the properties of the models in Table 2. For reference we list first the 
cooling-flow "CF" model: with no SMBH feedback (but allowance for SN feedback of types 
la and II, and taking into account the SMBH gravitational field), it allows the central SMBH 
to grow to mass far greater than observed in any system: clearly AGN feedback is needed. 



3.1. Purely radiative models 

We now discuss a couple of purely radiative models (Table 2, models RBO and RBO02), 
i.e. models in which the mechanical feedback effects are computed (following the Type B 
description) but not added to the hydrodynamical equations, and whose maximum radiative 
efficiency is e = 0.2 (RBO02) or e = 0.1 (RBO), values which bracket current observational 
estimates for the likely bolometric radiative efficiency (e.g., Fig. 2 in Yu & Tremaine 2002). 
We begin by noticing that we do not expect that these models will be as satisfactory as the 
models presented in CO07, at least for what concerns the final SMBH mass. In fact, we see 
from the first line in Table 1 that not only the dark matter halo is much more important at 
large radii than in our previous works (so that global degassing is more difficult), but also 
that in the central regions of the galaxy the gas is more tigthly bound. This is because the 
stellar density profile is now proportional to r~ 2 , so that the gravitational field is stronger 
in the central regions; in addition, the mass return from evolving stars is correspondingly 
more concentrated and the radiative losses more important. Finally, radiative feedback is 
reduced at low accretion luminosities, due to the ADAF prescription. This reduction has 
significant effects: in fact, it is during the low-luminosity phases that the stellar mass losses 
build up the galactic gaseous atmosphere responsible for the successive cooling catastrophe 
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episodes. A very low accretion luminosity during the quiescent phases contributes less to the 
global energetic balance, the subsoic outflow at large radii (mainly due to SNIa heating) is 
less favoured, and the time interval between the catastrophes is reduced (see also Fig. 8 in 
CO01). 

A first important result of the new models is that overall the main properties of the 
CO01 and CO07 models are confirmed, and episodic outbursts reaching I ~ 0.1 are common. 
After a first evolutionary phase in which a galactic wind is sustained by the combined heating 
of SNIa and thermalization of stellar velocity dispersion, the central "cooling catastrophe" 
of the galaxy gaseous halo commences, with the formation of a collapsing cold shell at ~ 1 
kpc from the center. In absence of the feedback from the central SMBH a "mini-inflow" 
would then be established, with the flow stagnation radius (i.e., the radius at which the 
flow velocity is zero) of the order of a few hundred pc to a few kpc: these decoupled flows 
are a specific feature of cuspy galaxy models with moderate SNIa heating (Pellegrini & 
Ciotti 1998). After the cooling catastrophe, the radiative feedback caused radiation pressure 
and radiative heating strongly affects the subsequent evolution, as can be seen in Fig. [2] 
where we show the luminosity evolution of the central AGN with time-sampling of 10 5 yrs. 
The bolometric luminosity (top panel) ranges between 0.001-0.1 of the Eddington limit (the 
almost horizontal solid line) at peaks in the SMBH output but, since obscuration is often 
significant, the optical accretion luminosity as seen from infinity can be much lower (bottom 
panel): we note that obscuration is commonly considered the explanation for red quasars 
(e.g., see Wang 2008). However, the central quasar is not always obscured and we see, in 
the lower panels of Fig. [21 that the optical luminosity exceedes ~ 10 44 erg s -1 in numerous 
bursts. As already found in CO07, the major AGN outbursts are separated by increasing 
intervals of time (set by the cooling time and by the secular decrease of the mass return rate 
from the evolving stellar population), and present a characteristic temporal substructure, 
whose origin is due to the cooperating effect of direct and reflected shock waves (from the 
inner rim of the spherical strongly perturbed zone): these outflowing shocks would be a 
likely place to produce emission of synchrotron radiation and electron and ionic cosmic rays, 
that are considered an additional source of feedback (Jiang, Ostriker & Ciotti, Paper V, in 
preparation; see also Sijacki et al. 2008). 

In model RB0 2 (left panels), at t ~ 6 Gyrs the SNIa heating, also sustained by a last 
strong AGN burst, becomes dominant and after a last major burst a global galactic wind 
takes place and the nuclear accretion switches to the optically thin regime. The SMBH mass 
accretion rate strongly oscillates as a consequence of radiative feedback, with peaks of the 
order of 10 or (more) M /yr, while during the final, hot-accretion phase the almost stationary 
accretion is ^10~ 2 M Q yr _1 , much less than the istantaneous mass return rate from the 
passively evolving stellar populations (M*<1M yr _1 ). In the optical, log(L|y t /LEdd) ~ 
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—5 (see Table 2), consistent with observations (e.g., see Ho 2008). However, the absolute 
accretion luminosity emitted (~ 3 10 43 erg s _1 ) is still very large when compared to some of 
the low-luminosity AGN observed (e.g., Pellegrini 2005a,b), showing that some additional 
form of feedback is still needed in the low-luminosity phase of the models. Note that in the 
last 6 Gyrs the SMBH virtually stops its growth, while the ISM mass first increases due to the 
high mass return rate of the evolving stellar population, and then decreases due to the global 
galactic wind induced by SNIa. During the entire model evolution, a mass ^3 1O 1O M of 
recycled gas has been added to the ISM from stellar mass losses. Approximately 2.1 x 1O 1O M 
has been expelled as a galactic wind, while ~ 1.4 x 10 10 M Q has been transformed into new 
stars, so that only 0.7% of the recycled gas has been accreted onto the central SMBH. We note 
that several observational indications exist supporting the idea that, while the majority of 
the stellar mass in elliptical galaxies may have formed at high redshifts, small but detectable 
star formation events (summing up to ^5 — 10% of the total stellar mass) may have occurred 
at low redshift (e.g., see Watabe, Kawakatu Sz Imanishi 2007; Pipino et al. 2008, Helmboldt, 
Walterbos & Goto 2008; Trager, Faber & Dressier 2008). Since the flow becomes nearly 
steady, the duty-cycle approaches unity for RB02, but for RB0 the short time in a strong 
burst at late epoch produces a small duty-cycle as observed in low redshift central SMBHs. 
Overall, the global behavior of these purely radiative models is very similar to the models in 
CO07, where we refer for details about the different evolutionary phases. Integrated values 
of their global quantities at the end of the simulation are listed in Table 2. Finally, it is 
important to stress that an identical model without SMBH feedback (i.e., e = in eq. [33]), 
but with the same star formation treatment of the model just described, produced a SMBH 
of (the much too high) final mass >1O 1O M , while the total mass in new stars was reduced to 
~ 3 x 1O 9 M (see model CF in Table 2). In addition, model CF does not present fluctuations 
in the starburst and ISM X-ray luminosities, thus showing the vital importance of SMBH 
feedback in the overall results. 

Of course, reducing the value of the radiative efficiency from 0.2 to 0.1, the number of 
bursts and the final mass of the SMBH increases. In particular, notice how major bursts 
persist to late times, (right panels in Fig. [2]): this behavior reduces considerably the derived 
value of /duty calculated at late times, and the various quantities of the model are reported 
in Table 2. 

Are these purely radiative models an adequate representation of observed AGN/elliptical 
galaxies? While far better than cooling flow models such as CF, which do not allow for 
radiative heating and radiation pressure caused by the AGN, they are clearly inadequate 
in several respects. At the end of the simulations, the ratio of the SMBH to stellar mass 
(also considering the added mass in new stars) is ~ 9 10 -3 (for model RB0) and ~ 4.7 10~ 3 
(for model RBO02), larger than the canonical ratio of 0.0013 (e.g., Yu & Tremaine 2002), 
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Fig. 2. — Purely radiative models. Dotted lines are the bolo metric accretion luminosity (top) 
and the optical SMBH luminosity corrected for absorption, i.e., as it would be observed from 
infinity (bottom), for the models RB0 2 (left panels), and RBO (right panels). Model RBO 
is identical to model RBO02, but with e = 0.1 instead of 0.2. As in CO07 at the center we 
fixed Lf Ujapt (Ri) = O.ILbh- In model RB0 02 quasar-like bursts at early times are followed 
at late times by a very quiet and passive accretion at late times with quite low optical 
luminosity. With a lower radiative efficiency gas is ejected less efficiently from the galaxy 
and intermittent bursts persist to late times. While both models are far better than models 
without radiative feedback in addressing the problems caused by standard cooling flows, 
neither acceptably match the known properties of AGNs/young stars in massive galaxies 
(see Table 2). 
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i.e. radiative heating alone does not sufficie to limit SMBH growth to the observationally 
allowed degree. The luminosity from young stars is ~ 10 40 ' 2 erg/s even at relatively low 
states between bursts, which would produce central regions which are far too blue compared 
to typical ellipticals and the AGN during these low states, while at only 10 _4 LEdd) is far 
more luminous (at ~ lO ia2 L ) than observed AGN in their low states (see, e.g., Pellegrini 
2005a,b, Wrobel, Terashima & Ho 2007; Gallo et al. 2007, Ho 2008). Of course, from Table 
2 it is apparent that a model with a peak radiative efficiency of eo = 0.1 is worst. 

Are these defects due to too low an assumed radiative efficiency? For example, note 
that in model RB0 with a peak efficiency of 0.1, the average (i.e., weighted over the mass 
accretion rate, see eq. [33]) EM efficiency of this model is seen to be 0.057 which is below 
that found by Yu & Tremaine (2002) by applying the Soltan (1982) argument. However, 
model RB0 2 (with a peak radiative efficiency of 0.2), has an average radiative efficiency 
of 0.124, but is still inadequate in representing real galaxies, especially at redshift zero. In 
particular the relatively high values of the ratio L^ u opt / L^dd combined with the final high 
values of Mbh produces central optical luminosities ^10 43 erg s -1 , considerably above what 
is observed from nearby elliptical in the "off" state. For these reasons we believe that the 
defects of the purely radiative models cannot be "fixed" by simply increasing eo; radiative 
feedback alone is insufficient (for massive galaxies) in limiting cooling flow induced star 
formation and central SMBH mass growth. 



3.2. Purely mechanical models 

From the previous discussion, it follows that mechanical energy input is also required. In 
order to study the effects of mechanical feedback we now switch off radiation feedback (while 
radiative cooling is obviously maintained active), and we explore purely mechanical models. 
In some sense, these models are complementary to models RB0 and RB0 2- In a first assay, 
we add mechanical feedback in the simple fashion normally adopted in several investigations, 
as a fixed numerical efficiency for the input of mechanical energy and for the amount of gas 
carried out by the wind. These purely mechanical Type A models are referred in Table 2 as 
MA. We then present the purely mechanical models of Type B (MB models in Table 2), i.e. 
models in which the wind efficiency from the circumnuclear regions depends on the accretion 
luminosity, declining with declining values of I, as is expected for radiatively driven winds. 
In both families, the model identification number in Table 2 increases at decreasing wind 
efficiency. In this paper we describe the general properties of the purely mechanical models, 
while in Paper II we focus on more specific quantitative aspects of purely mechanical Type 
B models. 
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Fig. 3. — Purely mechanical Type A models, i.e. models with a constant nuclear wind 
mechanical feedback efficiency Luminosity evolution of models MAO (e^ 1 = 5 1CT 3 , left 
panels), and MA3 (e^ = 5 10~ 5 , right panels). Dotted lines are the bolometric accretion 
luminosity (top) and the optical SMBH luminosity corrected for absorption, i.e., as it would 
be observed from infinity (bottom). The global properties of Type A models are given in 
Table 2. 
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Fig. 4. — Purely mechanical Type B models, i.e. models with a wind mechanical efficiency 
dependent on the (normalized) accretion luminosity I. Luminosity evolution of models MBO 
(the purely mechanical model associated with model RBO, left panels), and MB3 (right 
panels). The nuclear wind mechanical feedback efficiency at peak are = 5 10~ 3 and 
= 3 10~ 4 , respectively. Dotted lines are the bolometric accretion luminosity (top) and the 
optical SMBH luminosity corrected for absorption, i.e., as it would be observed from infinity 
(bottom). The global properties of Type B models are given in Table 2. 
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3.2.1. Type A models: constant mechanical efficiency 

In these models the constant wind mechanical efficiency ranges from = 5 10~ 3 in 
model MAO, as utilized for example by Di Matteo et al. (2005) and Johansson, Naab & 
Burkert (2008), down to = 3 10 -5 in model MA4. The resulting global properties of the 
whole family are reported in Table 2. 

Model MAO (Fig. [3j left panels), with the highest assumed steady efficiency is found 
in a state of permanent, global wind. The fact that the accretion luminosity is very low 
during the entire evolution is revealed by the very low averaged value of < €em >— 2.7 10~ 3 
obtained from eq. (33), a consequence of the ADAF implementation. With the assumed 
high mechanical efficiency we find that (consistent with other investigators) almost all the 
recycled gas produced by stellar evolution is ejected from the galaxy, and a very negligible 
amount of mass is added to the central SMBH. In addition, also star formation is maintained 
at very low levels. The final X-ray thermal luminosity of the model is far too low, orders of 
magnitude below that seen from normal giant elliptical galaxies (e.g., O'Sullivan et al. 2003). 

As we reduce e^f (models MAI and MA2), the evolution becomes more interesting. In 
fact, gas accumulates via the usual thermal instability, falls to the center and there generates 
an enormous explosion and, in a single burst, the galaxy is essentially evacuated of gas. 
There are no further bursts and the system is in permanent wind state until z = 0. This 
is essentially what is found by Di Matteo et al. (2005) with the only difference being that 
the gigantic outburst in their case is fueled by a gas rich merger while in our case it is fed 
by recycled gas deposited at the center via the classic Field thermal instability. Since all of 
these models (MAO, MAI, and MA2) have no bursts at late times and are "on" at a nearly 
constant but low level, the corresponding duty-cycles at late times are near unity; again, 
these models are too faint as thermal X-ray sources. 

In model MA3 (Fig. [3j right panels), the mechanical feedback is quite weak, and the 
galaxy presents a series of bursts, well separated by an increasing amount of time. The final 
SMBH mass is almost doubled, and also the new stellar mass is now significant, summing up 
to several 1O 9 M . The situation is even more extreme in model MA4. These models match 
observations in one respect. In fact, the occurrence of numerous bursts at late times is 
revealed by small duty-cycles (10~ 3 — 10" 4 ) so that they would be seen as AGN at late times 
with approximately that probability, and this is roughly consistent with observations. An 
inspection of Table 2 also reveals that AGN feedback can be important not only to suppress 
star formation when co-operating with SNIa to remove the gas from the galaxy, but also 
to induce star formation, especially in the central regions of the galaxies. This finding was 
already obtained in CO07, and proposed as a possible way, in alternative or in addition to 
galaxy merging (e.g., Hopkins et al. 2008a,b), to produce the central "extra-light" observed 
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in some elliptical galaxy (see Lauer et al. 2005). 

Clearly, these model cannot represent the cosmological evolution of galaxies, as all 
galaxies would be found systematically at very low values of their ISM X-ray luminosity, 
just consistent with discrete sources, but in net disagreement with observations if considering 
high efficiency models. In addition, a reduction of the fixed mechanical efficiency leads to 
a quite strong transition to the situation in which a significant amount of mass is added to 
the central SMBH. In other words, purely mechanical models with constant wind efficiency 
appear to be quite "stiff" in their behavior, switching from a strong global wind to a recurrent 
and significant bursting activity when reducing the efficiency. Extreme "fine tuning" would 
be required to produce an acceptable model and that result would probably not be robust 
to the normal statistical variations of parameters describing ellipticals in the Fundamental 
Plane. Due to the computational time required, such statistical investigation is however not 
attempted here, and it is postponed to Paper II. 

Guided by these results, we reduce the importance of feedback by considering a more 
plausible variation of the purely mechanical models (Type B models), in which the mechan- 
ical efficiency is a function of accretion luminosity, so that during bursts e w is moderate, and 
drops to very low values in the low-luminosity phases between bursts, in better accord with 
theoretical models for radiative accretion of AGN winds. 



3.2.2. Type B models: luminosity dependent mechanical efficiency 

We now move to describe the results for the family of MB models. As in the case 
of MA models, the value of the (peak) mechanical efficiency decreases from model MBO 
to model MB4. In particular, the peak values e^ 1 are coincident with the constant values 
of e w for the pairs MA0-MB0 and MA2-MB4, so that in these models the effects of the 
luminosity-dependent mechanical feedback efficiency can be studied by direct comparison. 

In general, the effect of the luminosity dependence of e w in Type B models is an increase 
in the final mass added to the central SMBH, as can be seen from Table 2. For example, at 
the end of the simulations AM BH ~ 1O 7 ' 17 M in model MAO, while AM BH = 1O 8 - 98 M in 
model MBO. Similarly, this is ~ 1O 7 75 M in model MA2 and 1O 9 - 44 M in model MB4. Note 
that these differences are mainly produced during the low-luminosity phases, as at the peaks 
the mechanical efficiency is identical to the corresponding MA models, and the comparison 
of the luminosity evolution in the right panels of Fig. [3] and Fig. H] illustrates this difference. 

As an additional example, in the right panels of Fig. H] one can see a large number 
of bursts during all the evolution of model MB3 (e^ 1 = 310~ 4 ), while model MAI (e^f = 
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2.5 10~ 4 ) just presents a single burst. And the accreted mass of the central SMBH (AM B h) 
is in far better agreement with observations (cf. MB2, MB3) than it is in the models with 
constant accretion efficiency. The present-day luminosity ratios of the (optical) accretion 
luminosity to the Eddington luminosity are much smaller than in the purely radiative models, 
but higher than in Type A models. The duty-cycles (Table 2, column 13) are in general quite 
small, due to the recurrent bursting activity which is still present in the models at late epochs. 
It is interesting that the duration of a single burst is of the order of 10 6 yrs, in accordance 
with observational estimates based on proximity effect (e.g., Kirkman & Tyler 2008); we also 
note how the luminosity averaged EM efficiencies, both in Type A and Type B models, are 
in the range 0.04< < e EM > ^0.08, in nice agreement with observational studies. 

The amount of new stars formed is of the same order as in Type A models, while the 
final mass of the ISM (and its X-ray luminosity) is larger, due to the less efficient feedback 
during the low luminosity phases. Also, the final SMBH masses are larger than those that 
are observed, and larger than in Type A models. Therefore, while in Type A models almost 
no mass is added to the central SMBH after the process of galaxy formation (and so no 
cosmological evolution in the Mbh — cr relation is expected in absence of merging, in Type B 
models a mass comparable to the initial SMBH mass is added on a cosmological time, and 
a consequent evolution of the M B h — cr relation is expected. We note that observations are 
now becoming able to probe such evolution (e.g., see Woo et al. 2008). Finally, during the 
low-luminosity phases, the nuclear emission is larger than what is observed in (very) low- 
luminosity AGNs, and this is due to the reduced contribution of the luminosity-dependent 
mechanical wind efficiency in Type B models. 

We then conclude that an additional form of feedback is needed both during the high 
luminosity phases (and this feedback is in the form of radiation, explored in Paper III, 
which is devoted to the study of combined models), and during the quiescent low-luminosity 
phases at late times. Obviously, this second form of feedback cannot be provided by radiation 
pressure and radiative heating as in the peak phases, and it is presumably provided by jets 
and/or thermally driven nuclear winds, not included in the present models. In sum, Type 
B models (ADAF-like) are more satisfactory than the normally computed (fixed efficiency) 
Type A models, but are still observationally inadequate. 



4. Discussion and conclusions 

In this paper we have investigated, with the aid of 1-D hydrodynamical simulations 
with unprecedented spatial and time resolution (where the cooling and heating functions 
include photoionization and Compton effects, and restricting to models which include only 
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Note. — Masses are in units of Solar Masses and luminosities in erg/s. In models of class A the wind efficiency is maintained constant, i.e., e w = e™. In models of class 'b the value of 
ejj? is reached when Lbh = 2Z;Edd- I n models with the subscript 02 the maximum radiative efficiency in eq. (8) is set to eo = 0.2 instead of 0.1 as in the other models. Mean efficiencies 
are calculated according to equation (34). AM* is the total amount of star formed during the model evolution, AM W is the total amount of ISM lost at 10i? c and M gas the instantaneous 
amount of gas inside 10R e . i B *Jj t is the fiducial SMBH luminosity in the optical as would be seen at infinity after absorption, with t = O.ILbh at the first grid point (see 

CO07 for details). The luminosity weighted wind solid opening angle < Af2 w > (in units of 4n) is calculated according to eq. J34), and the reported value refers to the final time of each 
simulation; in model MAO the value is not reported as Lbh < 0.1Z/Edd over all the evolution. The duty cycle is calculated for the effective optical accretion luminosity as would be seen 
from infinity. The first number correspond to a Universe age of 4 Gyr, the second to 13.7 Gyr, and the time span used for the computation is At = t/2. 

a Figures corresponding to quantities calculated but not added to the hydrodynamical equations. 

b The cooling flow model CF was stopped at 8 Gyrs. 
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radiative or only mechanical feedback in the form of nuclear winds), the effects of radiative 
and mechanical feedback on the gas flows in isolated elliptical galaxies. In particular, after 
having explored a few purely radiative feedback models, we explored two sequences of models 
in which the feedback is purely mechanical, and the only radiative effects allowed in the 
simulations are those due to gas cooling. The first sequence (Type A models) has a fixed 
mechanical efficiency as in most prior investigations, and the second (Type B models) allows 
the efficiency to increase with increasing accretion rate as indicated by the radiatively driven 
outflow models. It is difficult to know, either from fundamental physics calculations or from 
existing observations, what values of the mechanical energy efficiency to adopt, so we have 
tried a range of values from 5 10~ 3 to 10~ 5 to see which, if any, provides acceptable results, 
consistent with what is known about the properties of elliptical galaxies. 

The investigation is in line with that of our previous papers, and the general framework 
is maintained. However, the galaxy models and some of the relevant input physics have been 
substantially improved and extended (see Table 1). As before, the recycled gas, arising from 
the evolving stars in the inner several kpc of the galaxy (assumed a giant elliptical), necessar- 
ily drives a classical radiative instability and a collapse towards the center of metal rich gas. 
As a consequence, a star-burst occurs and the central SMBH is fed. We confirm that steady 
accretion on SMBHs is only possible at very low Eddington ratios, and consistently with the 
work of other authors we note that independent of whether the feedback is mechanical or 
radiative no steady flow appears to be possible for Eddington ratios above ~ 0.01, and that 
whenever the luminosity is significantly above this limit both the accretion and the output 
luminosity is in burst mode. The details of how much gas is accreted on the central SMBH 
vs. consumed in stars vs. ejected from the center by energy input from the starburst and 
AGN, are one of the main outcomes of these models. Relevant quantitative properties of the 
models are presented in Table 2, while the general results can be summarized as follows: 

1) Radiative heating and radiation pressure on the ISM by photons emitted by the cen- 
tral AGN and by the starburst, without any mechanical input, greatly reduces the "cooling 
flow catastrophe" problem, but leads to results that are still defective as compared to de- 
tailed observations of local elliptical galaxies, in that the central SMBH would be too bright 
and too massive, and the galaxy would be too blue, due to repeated bursts of central star 
formation. 

2) Adding mechanical energy from an AGN wind with fixed efficiency can address the 
problems posed by pure cooling flows but does not give a solution that in detail satisfies the 
observations. If the chosen efficiency is large, then a giant burst and an explosive degassing 
of the galaxy occurs (consistent with the results of Di Matteo et al. 2005, Johansson et al. 
2008, and other investigators). Only a small growth of the black hole occurs before the gas 
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content of the galaxy drops to levels below what is observed in real elliptical galaxies and the 
systems at redshift z = are computed to have thermal bremstrahlung X-ray luminosities 
orders of magnitude lower than those typically seen in nearby ellipticals. Finally, we do not 
get any late time AGN, and the duty cycle approaches unity. By contrast, very low values 
of the assumed efficiency do not prevent the cooling flows from causing more late time star 
formation (blue stellar cores) than is seen in most nearby ellipticals. It appears that there 
is no range of intermediate efficiencies that allows one to escape from both of these defects. 

3) Models with mechanical energy efficiency proportional to the luminosity, as indicated 
both by observations and detailed 2-D hydrodynamical simulations for radiatively driven 
winds, perform better, but are still inadequate. We thus conclude that mechanical energy 
input of the type adopted here - by itself - is unable to provide appropriate levels of feedback 
that would leave ellipticals at the current epoch with the properties that they are observed 
to have. In fact, we find clear evidence that two additional different forms of feedback are 
needed. The first, taking place during the high luminosity phases, it is in the form of radiative 
feedback, and is required to limit the growth of the central SMBH. The other form of missing 
feedback is required during the quiescent, low- luminosity accretion phases (in particular at 
late epochs), as revealed by accretion luminosity values that, altough very sub-Eddington, are 
still larger than observed. Of course, standard radiative feedback is not effective during such 
phases, and presumably the further reduction is provided by nuclear jets and/or thermally 
driven winds. We remark again that these results are based on ID simulations. It would 
be very important to have similar 2D or 3D simulations, with the same input physics, for 
a verification of the conclusions reported above, in particular when considering the issues 
related to stronlgy aspherical phenemona (such as the effective coupling of jets and nuclear 
winds with the ambient ISM), or the effect of the ISM angular momentum on the accretion 
rate. 

4) We stress that in our study the initial conditions represent a galaxy with a central 
SMBH closely following the Magorrian relation, and with a stellar population already formed: 
such a system would be called a "dead and red" galaxy, without much evolution expected. 
However, one of the general consequences of our exploration is the fact that the recycled 
gas from dying stars is a major source of fuel for the central SMBH, and this can induce 
substantial QSO activity, even in the absence of external phenomena such as galaxy merging, 
and the simulations reveal how complex the evolution of an isolated galaxy, subject to internal 
evolution only, can be. We note that recently similar conclusions have been reached, by using 
different considerations, also by other authors (e.g., see Pierce et al. 2007, Li et al. 2008, 
Kaufmann & Heckman 2008). 

Summarizing, neither models with purely radiative feedback nor those with mechanical 
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feedback as adopted here alone seem to be able to match all of the basic observations of 
elliptical galaxies (see also Sijacki et al. 2008). Models with purely radiative feedback are 
unable, for any reasonable radiative efficiency, to forestall sufficiently the collapse of gas onto 
the central black hole. While they do reduce the mass far below what it would have been in 
the no feedback ("cooling flow") case, they do not reproduce the observed Magorrian relation 
and leave the central black hole with masses greater than 1O 9 M , perhaps a factor of four too 
large. The models with purely mechanical feedback can, however, for an appropriate choice of 
the mechanical efficiency, match the well established M^—a relation, but only for efficiencies 
chosen to be so large as to leave the galaxies with much less gas than is actually seen in 
normal ellipticals as determined by their X-ray luminosities. A more detailed discussion of 
other aspects of the difficulties encountered by the models with solely mechanical feedback, 
the description of the most successful combined models (which utilize both radiative and 
mechanical feedback), and their observational properties, is kept for the following papers 
(Paper II, III, IV, and V). 

We thank Jenny Green, Jeremy Goodman, Silvia Pellegrini, Eve Ostriker and Jim Stone 
for useful discussions, and Annibale D'Ercole and Dennis McRitchie for helpful sugges- 
tions on coding issues. We also thank the anonymous Referee for a careful reading and 
for comments that improved the presentation of the paper. D.P. acknowledges support by 
the National Aeronautics and Space Administration under Grant/Cooperative Agreement 
No. NNX08AE57A issued by the Nevada NASA EPSCoR program. 

A. A physical basis for the phenomenological mechanical feedback equation: 
input from AGN wind/jet to hydrodynamic flow 

In this Appendix we describe the physics behind the phenomenological differential equa- 
tion (l2"5j) modeling the wind and jet mechanical interaction with the ISM of the galaxy: note 
that the treatment is formally identical for the jet and for the wind. 

How can we estimate the rate of energy deposition from a conical wind or jet emitted 
by a central AGN into the ISM of the galaxy in which the SMBH is embedded? We will 
take from Sects. 2.2 and 2.3 the estimates for the total outflowing mass rate in the wind 
or jet (M dw and M- v eqs. [T7] and [15] ; we denote these quantities as M wj ), and the total 
mechanical energy flux (Law and Lj, eqs. [20] and [21]; denoted as £ w j)- The solid angle of 
each cone, that we call in general AQ C (constant for the jet, and given by eq. [31] for the 
wind), has a linear opening angle C , so that AQ C = 7rG^. From elementary geometrical 
considerations, one obtains for the linear opening of the cone the relation R c (r) = r0 c (r), 
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M W j = t> w j(r)7ri?^(r)p W j(r), and L w j = 0.5M w jf^(r). As described in Sect. 2.3, the jets were 
taken to be very narrow and nearly relativistic, with energy efficiency increasing as Mbh 
decreases. The winds were taken to have much larger opening angles, much lower (but still 
supervirial) velocities and to have an energy efficiency that increased with increasing Mbh(^) 
(see Fig. 1 for details). 

Certainly the simplest description that one could propose would be to imagine that no 
energy or mass was transferred from the out-flowing supersonic flow to the ambient fluid 
until the internal momentum flux per unit area (declining as r -2 ) reached the level where it 
locally matched the pressure in the ambient fluid, and define that radius (-R w j) to be that of 
the working surface where the wind or jet dumped all of its mass, momentum and energy. 
In this simplest of all models there would be no input from the wind or jet until this point 
was reached as the jet/wind was driven loss-lessly through the ISM of the galaxy. So, below 
the working surface, where (3{r) = Pism/-Pwj < 1, the flow is conservative, but then as the 
pressure ratio approaches unity the mass and energy are deposited rapidly into a thin shell. 
While this approach has the virtue of simplicity and conservation of the vital quantities (and 
it was the first that we explored numerically), it is numerically unattractive and unstable. 

Next one could adopt the similar but smoother scheme described in eq. f l29|) but without 
the time dependent term. In this case there are some losses below and above the working 
surface, but the bulk of the mass and energy are dumped into the ambient medium over 
approximately one scale- length in radius in the vicinity of i? w j. This is the equation that 
we adopted in this paper for adding energy and mass from the wind/jet to the ambient 
flow. This approach has assumed implicity that all wind/jet outflow velocities are highly 
supersonic so that the quantities in the outflow at distance r from the center and at time t 
are essentially the same as at R\ and t. 

But it is straightforward to allow for propagation effects the finite time it takes for the 
wind/jet to flow from the first active grid point R\ to r. Since the total mass, momentum and 
energy flux (integrated over the sphere) are conserved beyond Ri aside from losses calculated 
explicitly (sinks, indicated in the following with S), we can write the more general equation 

dY dY 

- w + v^— = -S(r,t). (Al) 

We then finally choose to write the sink term as S(r,t) = Y/r(r,t), so that dt/r(r,t) is the 
fraction of the outflow mass (or energy or momentum) flux that is deposited in time dt at 
(r, t). Consistent with the concept of depositing the mass, momentum and energy smoothly 
in the vicinity of the working surface, we take 
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^(3(r,t), (A2) 
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where (3 is the pressure ratio. Inserting eq. (IA2I) in eq. flAlj) . and after some rearrangement, 
we obtain eq. (129]) . It reduces to our first, most simplified proposal (no explicit time de- 
pendent term) in the limits that either the flow is steady or it is of very large velocity (i.e., 
f w j max(t>, Cg)). The full equation (j2U|) . which allows for the lags in a time- dependent flow, 
provides more stable solutions but is more costly to implement than the non time-dependent 
form because the Courant condition for the time-step must be modified to be 

/ Ar Ar\ , 

At < mm (y?T^J' (A3) 

where Ar is the spatial grid spacing (for simplicity the contribution of artificial viscosity is not 
reported in equation above). In this paper we used the more economical time- independent 
version rather than the more accurate but costly to implement full version, which is used 
and discussed in Paper III. 

Now let us back up and ask if, for a steady conical flow, the time independent equation 
can be derived by any more formal arguments drawn from first principles. For a steady 
conical flow, which interacts with the ambient flow only at its boundary, it is easy to derive 
the following equations from the conservation laws. If the linear opening angle of the cone 
is C , then from geometrical considerations 

dln9 c _ c s 
d In r Oc^wj 

where the velocity of each fluid element in the outflow perpendicular to the radial direction 
is assumed to be the speed of sound c s in that flow at radius r. The equation of mass 
conservation is 

d\np = _ 2 _ d\\iv w] _ ^ d\n® c 

dlnr dlnr dhxr ' 

and we are now allowing the flow velocity to be a function of radius. The First Law of 
Thermodynamics (dU = dQ — PdV) can be written as 

2 dinc, = _ dlnp + j-l M {m) 
dlnr dlnr ct: dhxr 

where q(r, t) is the heat energy added per unit mass of the flow. The conservation of energy 
(neglecting gravitational energy in this super-virial flow) can be written as 



dlnr 7 — 1 v ^ d In r 



+ P(r)—^, (A7) 



where we have allowed for the work done on the external fluid, as described beelow. We have 
not required that the pressure in the wind/jet balances the external pressure and in such a 
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highly supersonic flow this would not naturally be the case (the jet fluid tending to rapidly 
cool due to the adiabatic expansion, i.e. c£(r)p(r) <C Pism)- We can define the equilibrium 
internal speed sound, c eq , that would correspond to pressure equilibrium across the cone 
boundary as 

p(r)c 2 cq (r) = P ISM (r), (A8) 

and in general we expect that, if c s (r) <C c cq (r), then the wind/jet will be unstable to 
rippling, internal shocks will develop, and we can plausibly assume that the heat generated 
by such processes will result in energy being taken from the bulk flow (reducing v w j) and 
being added to the thermal energy in the wind/jet, with the heat generated by the instability 
proportional to Cg q (r) — c 2 (r). Of course, additional phenomena such as strong shear and 
large density contrast at the wind/jet boundary are likely to give rise to instabilities there, 
regardless of the value of c s . Then, from dimensional analysis we would obtain 



— — = constant x — x \c'~-' { 
dhir u - 1 



constant x x [c 2 (r) - c 2 (r)]. (A9) 
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Note that the instabilities converting forward kinetic energy in the wind/jet to thermal 
energy, are also responsible for the source terms of mass, momentum and energy transferred 
from the wind/jet to the ISM along the boundaries of the cone. The set of six equations 
(A4)-(A9) would allow one to compute the radial dependence of the six variables p, f w j, 9 C , 
c eq , c s , q (or in principle their more complicated time- dependent versions) to determine the 
rate at which energy and mass are drained from the sides of the wind/jet and from its working 
surface. If however the instabilities are very efficient, we can short circuit this process, drop 
the q(r) variable and simply assume that heat is supplied to the outflow (and taken from the 
bulk motion) at the rate required to maintain c s (r) = c eq (r). Then, omitting several steps 
of the straightforward algebra, and approximating the flow by power- law solutions (i.e., all 
variables are taken to scale as powers of the radius r), one obtains, in the case of a steady 
wind/ jet 



d\nY 

-7j = -P( r ) 

dmr 



2 1 + /j 



i it \ ■ P= ~u ' ( Al °) 

7 — 1 + 7p(rj J dmr 

which, aside from the bracket is the same as the time-independent version of equation 
that we use in this paper. Typical values for the quantities in the brackets, considerably 
below the working surface are 7 = 5/3, \i = —2, (3 — > 0, so that the term in square brackets 
tends to the value of 2, giving a result which is very close to our time-independent version of 
eq. (29). The essential conclusion of this section is that, while a full 2D or 3D, time-dependent 
treatment of the energy and mss deposition from an out-flowing wind/jet is highly desirable 
and in fact essential for a careful, quantitative treatment of this problem, nevertheless, the 
simple formulation presented in eq. (29) does offer a conservative scheme which captures 
several of the basic features of the problem. 



-36- 



REFERENCES 

Antonuccio-Delogu, V. & Silk, J. 2008, MNRAS, 389, 1750 
Arav, N. et al. 2007, ApJ, 658, 829 

Begelman, M.C., & Nath, B.B. 2005, MNRAS, 361, 1387 

Begelman, M.C., & Ruszkowski, M. 2005, Phil.Trans. of Roy.Soc, part A, 363, n.1828, 655 
Bertin, G., et al. 1994, A&A, 292, 381 

Binney, J. 2001, in "Particles and Fields in Radio Galaxies Conference", ASP Conference 
Proceedings, Robert A. Laing and Katherine M. Blundell eds., vol. 250, p. 481 

Binney, J., & Evans, N.W. 2001, MNRAS, 327, L27 

Binney, J., & Tabor, G. 1995, MNRAS, 276, 663 

Blustin, A.J., Kriss, G.A., Holczer, T., Behar, E., Kaastra, J.S., Page, M.J., Kaspi, S., 
Branduardi-Raymont, C, & Steenbrugge, K.C. 2007, A&A, 466, 107 

Burkert, A., & Silk, J. 2001, ApJ, 554, L151 

Cappellari, M., et al. 2006, MNRAS, 366, 1126 

Castor, J.I., Abbott, D.C., & Klein, R.I. 1975, ApJ, 195, 157 

Cavaliere, A., & Vittorini, V. 2002, ApJ, 570, 114 

Chartas, C, Brandt, W.N., & Gallagher, S.C. 2003, ApJ, 595, 85 

Churazov, E., Sazonov, S., Sunyaev, R., Forman, W., Jones, C, & Bohringer, H. 2005, 
MNRAS, 363, L91 

Ciotti, L. 2009, La Rivista del Nuovo Cimento, 32, n.l., 1 

Ciotti, L., & Ostriker, J.P 1997, ApJ, 487, L105 (C097) 

Ciotti, L., & Ostriker, J.P. 2001, ApJ, 551, 131 (CO01) 

Ciotti, L., & Ostriker, J.P. 2007, ApJ, 665, 1038 (CO07) 

Ciotti L., Morganti, L., & de Zeeuw, PT. 2008, MNRAS, 393, 491 

Cowic, L.L., Ostriker, J.P, & Stark, A.A. 1978, ApJ, 226, 1041 



-37- 



Crenshaw, D.M., Kraemer, S.B., & George, I.M. 2003, ARAA, 41, 117 
Croton, D.J., et al. 2006, MNRAS, 365, 11 

Czoske, O., Barnabe', M., Koopmans, L.E.V., Treu, T., Bolton, A.S. 2008, ApJ, 384, 987 

de Zeeuw, P.T. 2001, in the Proceedings of the ESO workshop 1999 L. Kaper, E.P.J, van 
den Heuvel and P. A. Woudt Eds., pag. 78 

Diehl, S., & Statler, T.S. 2008 ApJ, 687, 986 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604 
Dorodnitsyn A., Kallman T., & Proga, D. 2008, ApJ, 675, 5 
Douglas, N.G., et al. 2007, ApJ, 664, 257 
Dubinski, J., & Carlberg, R.G. 1991, ApJ, 378, 496 

Dye, S., Evans, N.W., Belokurow, V., Warren, S.J., & Hewett, P. larXiv:0804.4002l 
Faber, S.M., et al. 1997, AJ, 114, 1771 
Fabian, A.C. 1999, MNRAS, 308, L39 

Fabian, A.C, Thomas, PA., Fall, S.M., & White III, R.E. 1986, MNRAS, 221, 1049 
Ferrarese, L., & Ford, H. 2005, Space Science Reviews, 116, 523 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Fukushige, T., & Makino, J. 1997, ApJ, 477, L9 

Gallo, E., Treu, T., Jacob, J., Woo, J.H., Marshall, P., & Antonucci, R. 2008, ApJ, 680, 154 
Ganguly R., & Brotherton, M.S. 2008, ApJ, 672, 102 

Gavazzi, R., Treu, T., Rhodes, J.D.. Koopmans, L.V.E., Bolton, A.S., Buries, S., Massey, 
R.J., & Moustakas, L.A. 2007, ApJ, 667, 176 

Gebhardt, K., et al. 2000, ApJ, 539, L13 

Ghisellini, G. larXiv:0807.2"640l 

Goncalves, T.S., Steidel, C.C., & Pettini, M. 2008, ApJ, 676, 816 

Graham, A.W., Erwin, P., Caon, N., & Trujillo, I. 2003, Rev.Mex.A.A., 17, 196 



-38- 

Granato, G.L., De Zotti, G., Silva, L., Bressan, A., & Danese, L. 2004, ApJ, 600, 580 
Haiman, Z., Ciotti, L. & Ostriker, J.P. 2004, ApJ, 606, 204 

Hamann, F., Kaplan, K.F., Rodriguez Hidalgo, P., Prochaska, J.X., & Herbert-Fort, S. 2008, 
MNRAS, 391, L39 

Helmboldt, J.F., Walterbos, R.A.M., & Goto, T. 2008, MNRAS, 387, 1537 
Hernquist, L. 1990, ApJ, 356, 359 
Ho, L.C. 2008 ARAA, 46, 475 

Holt, J., Tadhunter, C.N., & Morganti, R. 2008, MNRAS, 387, 639 
Hopkins, P.F., & Hernquist, L. larXiv:0809.3789l 

Hopkins, P.F., Hernquist, L., Cox, T.J., Robertson, B., Di Matteo, T., & Springel, V. 2006, 
ApJ, 639, 700 

Hopkins, P.F., Cox, T.J., Dutta, S.N., Hernquist, L, Kormendy, J, & Lauer, T.R. 2008a 
laj5Gv:0805.3533l 

Hopkins, P.F., Lauer, T.R., Cox, T.J., Hernquist, L., & Kormendy, J. 2008b larXiv:0806.2325l 

Humphrey, P.J., Buote, D.A., Gastaldello, F., Zappacosta, L., Bullock, J.S., Brighenti, F., 
& Mathews, W.G. 2006, ApJ, 646, 899 

Jaffe, W. 1983 MNRAS, 202, 995 

Jaffe, W., Ford, H.C., O'Connell, R.W., van den Bosch, F.C., & Ferrarese, L. 1994, AJ, 108, 
1567 

Johansson, PH., Naab, T., & Burkert, A. 2008, Astron.Nachr. 329, 956 

Jolley, E, & Kuncic, Z. larXiv:0802.0"902l 

Kauffman, G., & Heckman, T.M. larXiv:0812. 12241 

King, I.R. 1972, ApJ, 174, L123 

King, A.R. 2003, ApJ, 596, L27 

Kirkman, D., & Tytler, D. 2008, MNRAS, 391, 1457 



-39- 



Koopmans, L.V.E., Treu, T., Bolton, A.S., Buries, S., & Moustakas, L.A. 2006, ApJ, 649, 
599 

KSnigl A. 2006, MemSAIt, 77, 598 
Kormendy, J., & Richstone, D. 1995, ARAA, 33, 581 
Kurosawa, R. & Proga, D. 2008a, ApJ, 674, 97 
Kurosawa, & Proga, D. 2008b, larXiv:0812.3T53l 
Lauer, T.R., et al. 2005, AJ, 129, 2138 

Li, C, Kauffmann, G., Heckman, T.M., White, S.D.M., & Jing, Y.P 2008, MNRAS, 385 
1915 

Magorrian, J., et al. 1998, AJ, 115, 2285 
Marconi, A. & Hunt, L.K. 2003, ApJ, 589, L21 
Martinez-Sansigre, A., & Taylor, A.M. arXiv:0810.3"920l 

Martini, P. 2004, in Coevolution of Black Holes and Galaxies, ed. L.C. Ho, 169 
McCarthy, I.G., Babul, A., Bower, R.G., & Balogh, M.L. 2008, MNRAS, 386, 1309 
McLure, R.J., & Dunlop, J.S. 2002, MNRAS, 331, 795 
McNamara, B.R, & Nulsen, P.E.J. 2007, ARAA, 45, 117 
Merloni, A. & Heinz, S. 2008 MNRAS, 388, 1011 

Milosavljevic, M., Bromm, V., Couch, S.M., & Oh, S.P arXiv0809.2404 
Morgan, C.W., Kochanek, C.S., Morgan, N.D., & Falco, E.M. arXiv0707.0305 
Murray, N., Quataert, E., & Thompson, T.A. 2005, ApJ, 618, 569 
Narayan, R. & Yi, I. 1994 ApJ, 428, L13 

Navarro, J.F., Frenk, C.S., & White, S.D.M. 1997, ApJ, 490, 493 
Nayakshin, S., & Sunyaev, R. 2005, MNRAS, 364, L23 

Nayakshin, S., Dehnen, W., Cuadra, J., & Genzel, R. 2006, MNRAS, 366, 1410 
Noble, S.C., Krolik, J.H., & Hawley, J.F. larXiv:0808.3140l 



-40- 

Omma, H., Binney, J., Bryan, G., & Slyz, A. 2004, MNRAS, 348, 1105 

O'Sullivan, E., Ponman, T.J., & Collins, R.S. 2003, MNRAS, 340, 1375 

Ostriker, J.R, & Ciotti, L. 2005, Phil.Trans. of Roy.Soc, part A, 363, n.1828, 667 (OC05) 

Parriott, J.R., & Bregman, J.N. 2008, ApJ, 681, 1215 

Pellegrini, S. 2005a, ApJ, 624, 155 

Pellegrini, S. 2005b, MNRAS, 364, 169 

Pellegrini, S., & Ciotti, L. 1998, A&A, 333, 433 

Pellegrini, S., Ciotti, L. & Ostriker, J. P. 2009, Advances in Space Research, in press 

Peterson, J.R., & Fabian, A.C. 2006, Phys.Rep., 427, 1 

Pierce, CM., et al. 2007, ApJ, 660, L19 

Pipino, A., Silk, J., & Matteucci, F. larXiv:0810.2045l 

Pipino, A., Kaviraj, S., Bildfell, C, Hoekstra, H., Babul, A., & Silk, J. larXiv:0807.2760l 
Prochaska, J.X., & Hennawi, J.F. larXiv:0806.0862l 
Proga, D. 1999, MNRAS, 304, 938 
Proga, D. 2007a, ApJ, 661, 693 

Proga, D. 2007b, in The Central Engine of Active Galactic Nuclei, eds. L.C. Ho, & J.M. 
Wang, ASP Conf. Ser., 373, 267 

Proga, D., & Kallman T. 2002, ApJ, 565, 455 

Proga, D., & Kallman T. 2004, ApJ, 616, 688 

Proga, D., Stone J.M., & Drew J.E. 1998, MNRAS, 295, 595 

Proga, D., Stone J.M., & Kallman T. 2000, ApJ, 543, 686 

Proga, D., Ostriker, J.P, & Kurosawa, R. 2008, ApJ, 676, 101 

Rafferty, D.A., McNamara, B.R., & Nulsen, P.E.J. 2008, ApJ, 687, 899 

Renzini, A., Ciotti, L., D'Ercole, A. & Pellegrini, S. 1993, ApJ, 419, 52 



-41 - 



Riciputi, A. Lanzoni, B., Bonoli, S., & Ciotti, L. 2005, A&A, 443, 133 
Rusin, D. & Kochanek, C.S. 2005, ApJ, 623, 666 
Rusin, D., Kochanek, C.S., & Keeton, C.R. 2003, ApJ, 595, 29 
Saglia, R.P., et al. 1993, A&A, 403, 567 

Sazonov, S.Yu., Ostriker, J.R, & Sunyaev, R. 2004, MNRAS, 347, 144 

Sazonov, S.Yu., Ostriker, J.R, Ciotti, L., & Sunyaev, R.A. 2005, MNRAS, 358, 168 

Sazonov, S.Yu., Revnivtsev, M., Krivonos, R., Churazov, E., & Sunyaev, R.A. 2007, A&A, 
462, 57 

Sazonov, S.Yu., Krivonos, R., Revnivtsev, M., Churazov, E., & Sunyaev, R.A. 2008, A&A, 
482, 517 

Schawinski, K., Lintott, C.J., Thomas, D., Kaviraj, S., Viti, S., Silk, J., Maraston, C, Sarzi, 
M., Yi, S.K., Joo, S.J., Daddi, E., Bayet, E., Bell, T., & Zuntz, J. larXiv:0809. 10961 

Shakura N.I., Sunyaev R.A. 1973, A&A, 24, 337 

Shi J. & Krolik J. 2008, ApJ, 679, 1018 

Sijacki, D., Springel, V., Di Matteo, T., & Hernquist, L. 2007, MNRAS, 380, 877 
Sijacki, D., Pfrommer, C, Springel, V., & Ensslin, T.A. 2008, MNRAS, 387, 1403 
Silk, J., & Rees, M.J. 1998, A&A, 331, LI 
Soltan, A. 1982, MNRAS, 200, 115 

Somerville, R.S. 2008, in Panoramic views of galaxy formation and evolution, ASP Confer- 
ence Series, vol.399, eds. T. Kodama, T. Yamada, K. Aoki, p. 391 

Springel, V., Di Matteo, T. & Hernquist, L. 2005, MNRAS, 361, 776 

Tabor, G.,& Binney, J. 1993 MNRAS, 263, 323 

Trager, S.C., Faber, S.M., Dressier, A. 2008, MNRAS, 386, 715 

Treu, T. & Koopmans, L.V.E. 2002, ApJ, 575, 87 

Treu, T. & Koopmans, L.V.E. 2004, ApJ, 611, 739 



-42- 



Turnshek, D.A. 1988, in QSO Absorption Lines: Probing the Universe, ed. J.C. Blades, D.A. 
Turnshek, & C.A. Norman (Cambridge: Cambridge Univ. Press), 17 

Vernaleo, J.C. & Reynolds, C.S. 2006, ApJ, 645, 83 

Voit, CM., & Donahue, M. 2005, ApJ, 634, 955 

Wang, J.M. 2008, ApJ, 682, L81 

Watabe Y., Kawakatu, N., & Imanishi M. 2008, ApJ, 677, 895 

Woo, J.H., Treu, T., Malkan, M.A., & Blandford, R.D. larXiv:0804.0235l 

Wrobel, J.M., Terashima, Y., & Ho, L.C. 2008, ApJ, 675, 1041 

Wyithe, J.S.B., & Loeb, A. 2003, ApJ, 595, 614 

Yu, Q., & Tremaine, S. 2002, MNRAS, 335, 965 



This preprint was prepared with the AAS IATgX macros v5.2. 



